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ABSTRACT 

Establishing the relative role of internally and externally driven mechanisms responsi¬ 
ble for disc and bulge growth is essential to understand the evolution of disc galaxies. 
In this context, we have studied the physical properties of disc galaxies without clas¬ 
sical bulges in comparison to those with classical bulges since 2 0.9. Using images 

from the Hubble Space Telescope and Sloan Digital Sky Survey, we have computed both 
parametric and non-parametric measures, and examined the evolution in size, concen¬ 
tration, stellar mass, effective stellar mass density and asymmetry. We find that both 
disc galaxies with and without classical bulges have gained more than 50% of their 
present stellar mass over the last ~8 Gyrs. Also, the increase in disc size is found to be 
peripheral. While the average total (Petrosian) radius almost doubles from z ~ 0.9 to 
z ~ 0, the average effective radius undergoes a marginal increase in comparison. Ad¬ 
ditionally, increase in the density of the inner region is evident through the evolution 
of both concentration and effective stellar mass density. We find that the asymmetry 
index falls from higher to lower redshifts, but this is more pronounced for the bulgeless 
disc sample. Also, asymmetry correlates with the global effective radius, and concen¬ 
tration correlates with the global Sersic index, but better so for higher redshifts only. 
The substantial increase in mass and size indicates that accretion of external material 
has been a dominant mode of galaxy growth, where the circumgalactic environment 
plays a significant role. 

Key -words: galaxies: bulges - galaxies: evolution - galaxies: high-redshift - galaxies: 
structure. 


1 INTRODUCTION 

One of the major challenges linked to the morphologi¬ 
cal evolution of disc galaxies is to understand the forma¬ 
tion and evolutioir of their bulge component. The bulge 
of a disc galaxy, observationally, is the central component 
which contains all the light that is in excess of an inward 
extrapolation of a constant scale-length exponential disk 
llWvse. Gilmore Franrj|l997l : lButal[201^ . Internal and ex¬ 
ternal mechanisms can explain the growth of bulges iir disc 
galaxies. 


Externally driven bulge growth in discs can oc¬ 
cur through major mergers, multiple minor merg¬ 
ers, as well as accretion of small comp onen ts or 


satellites JKauffmann. White &: Guiderdonil 19931: 

Baugh. Cole & FrenkI Il996l: 

Aguerri. Balcells & Peletied 

2OOII: Bournaud. Jog & Combes 

2OO7I: Hookins et al. 201(J). 


The fraction of baryons which lose angular-momentum 
due to the dynamical friction produced in these mergers 
and accretion processes fall to the centre of the galaxy 


to form the bulge component (IParry, Eke fc FrenkI l2009l : 
lGovernato_et^_ajjl20]ii_ Broo]£_et aljj201ll. see also recent 
review by Brooks &: ChristensenI 2015l l . 


In terms of internally driven, there can be two bulge¬ 
building mechanisms. The first one is known as secu¬ 
lar evolution. In this mechanism, slow, internally cre¬ 
ated processes contribute to the rearrangement of an¬ 
gular momentum and mass inside disc galaxies, lead¬ 
ing to the growth of pseudo-bulges, which hav e proper¬ 
ties s imilar to discs iTBinnev fc Tremaine! 198^ Combed 


ties s iinUar to discs llijinnev &: iremamel llbtoil: lUombesI 

200^ Komieirdy fcj^nnicuttll2004l : iKormendv et al.ll2010l : 


Saha fc GerhardI 2013l l. In the second mechanism, the co¬ 


alescence of giant star forming clumps, due to internal 

J ravitational instabilities, leads to bulge f ormation in discs 
Elmegreen. Bournaud fc ElmegreerJl2008| . see recent review 
by Bournaud 2015l l. The relative importance of these mech¬ 
anisms in forming present day discs with varied bulge prop- 
erties is as yet n ot quantitatively known (see discussion in 

lKormen^lAl5l ~). 
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One of the most prominent nnresolved issnes associ¬ 
ated with bulge evolntion is the conspicnous presence of 
a large number of bulgeless galaxies (disc galaxies with¬ 
out classical bulges) in the local universe. This is be¬ 
cause, naively, a relatively low number of disc galax¬ 
ies, formed as per the ACDM structure formation sce¬ 
nario, are expected to survive without forming a classi¬ 
cal bulge in the centre. The intense amount of merger vi¬ 
olence associated with hierarchical clustering is expected 
to destroy the fragile thin disc of stars. Thus, their pres- 


tions (iKautsch et al.ll2006l: iKautsch 

I 2 OO 9 I: iKormendt 

f et al. 

I 2 OIOI). models dHookins et al.l 20091: 

Governato et al. 

2010|), 

simulations |Scannapieco et al.l l2009l: Zavala et al. 

2012) 

and detailed in reviews (iBaushI l2006l: iBenson 

201f: 


Observational evidence is required to pin down the rel¬ 
ative role of various mechanisms in bulge growth and under¬ 
stand the formation of discs without bulges. The constraint 
involved is that for distant galaxies only the very basics of 
structure and morphology can be investigated. However, re¬ 
cent developments in observational facilities and in para- 
metri c and non-para metric measurements of galaxy struc¬ 
ture (IConselicd l2014l . and references therein) have enabled 
us to use galaxy structure to measure fundamental proper¬ 
ties of distant galaxies. These properties can then be com¬ 
pared with those of nearby galaxies to determine bulge-disc 
evolution. 

_ There has been some progress in this regard. iGadottil 

ll2009ll found an overlap in the structural properties of 
pseudo- and classical bulges indicating that the different pro- 
cesses for bulge grow t h mig ht have happened concomitantly. 
IParrv. Eke fc FrenkI (l2009l i studied the models based on 
ACDM cosmology and found that most spiral bulges acquire 
thei r stellar mass through minor merger s and disc instabili¬ 
ties. iMo. van den Bosch fc Whitel (120101 1 explained the for¬ 
mation of bulges through mergers, secular inte rnal processes 
as wel l as misaligned/perturbed infalling gas. IWatson et all 
(1201111 studied the neutral hydrogen properties of 20 bul¬ 
geless galaxies to compare the role of mergers vs internal 
processes in bulge formation. They report that even though 
some of the discs have distinct outer compone nts indicating 
recent interaction, the discs remain bulgeless. I Zavala et al.l 
(I 2 OI 2 II examined the impact of mergers on the growth of 
bulges using simulated data. They found that the main chan¬ 
nels of bulge mass assembly are stars from infalling satel¬ 
lites, and stars transfer red from primary di scs due to merger 
induced perturbations. IPerez et al.l (l2013l l showed through 
simulations that strong disc instabilities at high redshifts 
lead to classical bulge formation, which cannot be prevented 
by eve n the most energetic supernova feedback. iBruce et al.l 
(|2014| found that from redshift 3 to 1, galaxies move from 
disc dominated to increasingly bulge dominated morphol¬ 
ogy- 


Although these works have given us considerable insight 
into disc structural evolution, the relative role of internal and 
external mechanisms in the formation and evolution of discs 
of varied bulge types is as yet not established. To achieve 
this, it will be insightful to examine the evolution in the in¬ 
ner region properties of disc galaxies of different bulge types 
in a relative manner. The crucial aspect in that direction is 


that the separation of disc galaxies according to their bulge 
type has to be achieved in a quantitative and robust manner. 

Thus, in this paper, we undertake a comparative study 
of the evolution of disc galaxies with and without classi¬ 
cal bulges. Disc galaxies without classical bulges are, by 
definition, labelled as “bulgeless” and those with classical 
bulges are labelled as “normal”. The two morphological 
types are separa ted using both Sersic index and Kormendy 
relation c r iteria (lGadottill2009li . This work is a follow-up to 
ISachde\^ (l2013h . such that the galaxy sample is the same 
and we also utilize the parametric measures (through Sersic 
function fitting) derived in that paper. Here we derive the 
non-parametric measures (Petrosian radius, concentration, 
asymmetry) along with rest-frame colours, total stellar mass 
and effective stellar mass density. The evolution in the inner 
region properties is thus examined through both parametric 
as well as non-parametric measures. 

The study is done since 2 ~ 1, when the galax- 
ies have Just formed a familiar Hubble seq uence structure 
l|Gonselice et al.ll201ll : [Mortlock et al.ll2013l i. to the present 
epoch, where they have develop ed and settled i nto few dis¬ 
tinctly i dentifiable categories dGadottil I2OO9I : iButal I2OI3I : 
lGraharnll2013tl . This time interval has the potential to reveal 
the major processes involved in bulge and disc evolution. 

For distant galaxies we make use of deep imaging from 
the Great Observatories Origins Deep Survey (GOODS) 
obtained using Hubble Space Telescope (HSTjA Advanced 
Ga mera for Surveys in the Chandra Deep Field South (CDF- 
S) (iGiavalisco et al.1l2004l i. For local galaxies, images are 
from the NASA Sloan Atl afl, based on the Sloan Digital 
Sky Survey dBlanton et al.llioOSbIi . We obtain the paramet¬ 
ric and non-parametric measures for bright (Mb ^-20) disc 
dominated (Sersic index n < 2.5) galaxies in three redshift 
ranges (0.77^ z <1.0, 0.4< 2 <0.77 and 0.02< 2 <0.05) for 
rest-frame H-band. We then examine the evolution of the 
inner-region properties of the bulgeless disc galaxies in com¬ 
parison to the normal disc galaxies over the three redshift 
ranges. 

We consider a flat A-dominated universe with Ha = 
0.73, Hm = 0.27, Ho = 71 km sec“^ Mpc“^. In Section 2, we 
describe our data in terms of sample selection, preparation of 
images, defining and computation of various parameters. In 
Section 3, we present the results obtained by examining the 
evolution of size, concentration, stellar density and asym¬ 
metry. We also report the correlation seen between various 
parameters. In Section 4, we list the primary results of this 
work and discuss their implications in the light of previous 
studies. 


^ Based on observations obtained with the NASA/ESA Hubble 
Space Telescope, which is operated by the Association of Uni¬ 
versities for Research in Astronomy, Inc.(AURA) under NASA 
contract NAS 5-26555. 

^ http://www.nsatlas.org. Funding for the NASA-Sloan Atlas 
has been provided by the NASA Astrophysics Data Analysis Pro¬ 
gram (08-ADP08-0072) and the NSF (AST-1211644). 
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2 DATA 

2.1 Sersic parameters 


In a previous work l|Saclidevall201.ll 'l . images taken from HST 
- ACS in V (F606W), i (F775W), and z (F850LP) filters 


were used to obtain the rest-frame B-band properties of the 
galaxies lying in the CD F-S with redshift rangin g from 0.4 
to 1.0. First, SExtractor llBertin fc Arnout^llOOd i was used 
to iden tify the sou rces in the z-band image. Then, single 
Sersic (I Sersidll 968ll components were fit (using GALFIT, 
IPeng et al.l 2002l ~) on all the three filter images. The Sersic 
profile for the variation of a galaxy’s surface brightness from 
its centre is given as: 


7(r) = /e exp[-6„(( —- 1)], 

Te 


( 1 ) 


where n (the Sersic index) controls the degree of curvature 
of the profile, le is the surface brightness at r^, and bn is 
a constant such that is the half-light radius for a given 
value of n. We thus obtained parameters such as apparent 
total magnitude, half-light radius and Sersic index for all the 
galaxies in the three filters. 

Redshifts we re obtained from the COMBO-17 survey 
llWolf et al.ll2004li . The V filter provides rest-frame R-band 
properties for galaxies lying in the redshift range of 0.4 to 
0.6. The z filter provides the same for the redshift range of 
0.8 to 1.0. For the redshift range of 0.6 to 0.8, both i and 
2 filters can be used and thus an average of the properties 
obtained using both the filters was employed. The param¬ 
eters were thus obtained for 4124 sources in the rest-frame 
R-band for the full redshift range (0.4 to 1.0). 

Using the redshifts and accepted cosmological parame¬ 
ters, we then computed absolute magnitude, half-light radius 
in kpc and surface brightness in mag/arcsec^. The absolute- 
magnitude, M, for the galaxies is calculated using the rela¬ 
tion: 


M = m — 5 * logio{DL * 10^) — K, 


( 2 ) 


where Dl is the luminosity-distance in Mpc and K is the 
K-correction term that accounts for the difference between 
the observed band and rest-frame ba nd. It depends on the 
object’s spectral energy distribution (lOke fc Sandagel[l968l : 
iHogg et '^l2002h and, for a power-law continuum, it is given 
by the relation: 


Kcont = -2.5 * (1 -I- a„) * log(l -I- z), 


( 3 ) 


where is t he slope of the contin uum and has a canonical 
value of -0.5 ([Richards et al.|[2006fl . Therefore: 

Mb = m — 5 * log(R>i, * 10®) -I- 2.5 * log(v^l -|- z)). (4) 


The half-light radius of the galaxies that we obtain from 
Sersic component fitting is in pixels. They were converted 
into arc-seconds according to the plate-scale of the telescope 
and then into radians. The intrinsic half-light radius in kpc 
was then calculated using the relation: 


Re = Da* 1000 * Ae, (5) 

where Da is the angular-diameter-distance in Mpc and A0 
is the radians covered on the detector by the half-light ra¬ 
dius. 

The redshift-magnitude distribution of the galaxies was 
examined in the 0.4-1.0 redshift range for galaxies with 


Ms>-20 and Mb^-20 (shown in Fig. 1 of ISachdevall20l3 b 
For Mb>-20, galaxies with lower luminosities were not seen 
at high redshift ranges at all. However, for Ms ^-20, the 
number of galaxies was seen to be evenly distributed. Also, 
the number of galaxies in two equal comoving volume red¬ 
shift bins (0.4-0.77 and 0.77-1.0) was found to be almost 
the same for Ms^-20. Additionally, the magnitude limit for 
a reliable redshift estimate from COMBO-17 is mz~23.5, 
which for our upper redshift limit of 2 = 1.0 corresponds to 
Ms~-20. Based on the depth of HST imaging, and the red¬ 
shift accuracy limit of COMBO-17, a magnitude cut of -20 
was applied on the sample. We, thus, obtained 727 sources 
in the rest-frame R-band (0.4^ z <1.0) with Ms^-20 (elab¬ 
orated in ISachdevall20l3 ~). 


2.1.1 Separating bulgeless disc and normal disc galaxies 

The Sersic index value of 2.5 is employed in numerous stud¬ 
ies to separa te early-type (n>2.5) and late-type (n<2.5) 
galaxies (e.g . Ravindranath et al.ll2004l : iBarden et al.l[2005l : 
Ivan der Wefcoosi '). We used these criteria to obtain 496 late- 
type (or disc dominated) (n<2.5), bright (Ms ^-20) galaxies 
in the rest-frame R-band. 

The bulgeless disc galaxies include disc galaxies with¬ 
out bulges and those with pseudo-bulges. Pseudo-bulges are 
the bulges which have a higher ratio of ordered-motion 
to random-motion. Since they exhibit nearly exponential 
brightness profiles, disc galaxies with pseudo-bulges are, 
there f ore, considered as bulg eless llKormendv &: KennicuttI 
I 2 OO 4 I : iKormendv et al] l201oh . Thus, to separate bulgeless 
disc and normal disc galaxies, we separated discs with no¬ 
bulge or pseudo-bulge from discs with classical bulge in our 
disc dominated sample. 

The Sersic index values ranging from 1.7-2.0 have 
been suggested by many studies for t he separation of 


classical bulges from p s eudo bulges (e.g. Shen et all l2003l : 


lUaurikainen et al]l2007l : iFisher fc Drorvll2nn8l i. To obtain a 
Sersic index limit for our sample, we divided the entire sam¬ 
ple (i.e. without the magnitude cut, 4124 sources) into three 
ranges (0.8> n, 0.8^ n <1.7, 1.7^ n) with each range get¬ 
ting almost equal number of sources. We then examined the 
distribution of the mean half-light radius against the abso¬ 
lute ma gnitude bins for different Sersic index ranges (see 
Fig. 2 of Sachdevjl2013fl . 

We found that a value of n~1.7 divides galaxies into 
two groups where each group follows a particular half-light 
radius-magn itude (R-M) re lation, independent of n (shown 
in Fig. 2 of Sachd e^ 1201,31 ! . This is in striking agreement 


in fi g. Z ot lt>acndevall2Ul,-il l. inis is in striking agreement 
with I Shen et al.l (1200311 . who found the same value based 
on SDSS data, and claimed that the cut separates galaxies 
with exponential surface brightness profile (Sb/Sc) from the 
galaxies which do not have such profiles. 

In addition to applying the Sersic ind ex limit of n~1.7 , 
we have applied the Kormendy relation (iKormen d^Iilzl) 
to ascertain the separation of bulgeless disc and normal disc 
galaxies. This is based on the fact that, since the Kormendy 
relation is followed by elliptical galaxies, galaxies which are 
bul geless should s how themselves as outliers to the relation 
(see lGadottill2009l . for details). 

To achieve this, we fitted a linear relation to the sur¬ 
face brightness vs log-size data of elliptical galaxies in our 
sample, i.e. those with 2.8 ^ n < 4.5 from the 727 sources 
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(0.4^ z <1.0, Mfl^-20). Those disc galaxies (n < 2.5) which 
were lying below the ±3cr value of the zero-point (with fixed 
slope) were taken as outliers. The outliers obey the following 
relations: 

Me.B > 19.36 + 2.92 log(i?e,s), 0.4 ^ 2 : < 0.77 ( 6 ) 

and 

> 19.32 + 2.92 log(7?e,B), 0.77 ^ 2 < 1.0. (7) 

More than 80% of the disc galaxies which were found to 
be bulgeless according to these relations were also seen to 
have n<1.7. The two criteria are thus complementary to 
each other. Only those galaxies which satisfied both criteria, 
i.e. had Sersic index less than 1.7 and were outliers to the 
Kormendy rela tion, were chosen to be bulgeless (Fig. 3 of 
ISachde^^l2013^ . The bulge/total light ratio was found to be 
less than 0.2 (or 20%) for our bulgeless sample. This process 
of morphological determination was found to separate the 
galaxies in a similar mann er in the infrared as they do in 
the optical (ISachdevall2013l i. 


2.1.2 Overall sample obtained 


We obtained Sersic parameters in rest-frame B-band for 
496 bright, disc dominated galaxies separated into 186 bul¬ 
geless disc galaxies and 310 normal disc galaxies, in two 
equal comoving volume redshift ranges (0.77^ 2 : <1.0 and 
0.4^ 2 : <0.77). 

In addition to this, a low redshift (0.02^ 2 <0.05) 
catalog of dis c dominated galaxies was taken from the 
NYU-VAGC (jBlanton et al.l [2005bl l3l to establish a lo¬ 
cal sample. We obtained rest-frame B-band parameters 
from the single Sersic fit parameters in g and r fil- 
ters. Their AB magnitudes are galac tic-extinction corrected 
( Schleeel. Finkbeiner fc DavisI Il998l l and also K-corrected 
I Blanton et al ■II2003I) to the rest-fra me b andpasses. We used 
relations from Fukugita et al.l (|l996l l an d l.Tester et al.l (120051 1 
to obtain the absolute magnitude in rest-frame B-band: 


g' = V + 0.56 *{B-V)- 0.12, 


(8) 


r'= 1/- 0.49 *(B-F)-f 0.11. (9) 

Using those two equations we get: 

B = 1.419 * 5 '- 0.419 *r'-h 0.216. (10) 

The Sersic half-light radii of the g and r-bands were con¬ 
verted from arc-secon ds to kpc according to their redshift. 
The relati ons given Iw Barden et al.l (12005 1. that were found 
using the Ide Jond (Il996ll data, were used to obtain half¬ 
light-radii in the rest-frame B-band. Using the NYU-VAGC 
catalog, they studied the ratio of half-light sizes in the five 
SDSS bands to the size measured in one band as a function 
of wavelength. The results were accurately described by a 
linear fit (only ±3% correction factor ) and we r e in s triking 
agreement with the fit found for the Ide Jonj (Il996l l data. 
The slope gives average corrections to obtain the rest-frame 
sizes: 


Re{V) = 1.011 *Be(r), 


( 11 ) 


Be(B) = 1.017 *Be(U). 


( 12 ) 


Using those two we get: 

Be(B) = 1.017* 1.011 >i<Be(r). (13) 

We obtained the final local catalog of 764 galaxies in 
rest-frame B-band, in the redshift-range 0.02^z<0.05 with 
Mb ^-20. After that we applied the Sersic index criteria to 
separate disc galaxies. Then disc galaxies were further sepa¬ 
rated into a sample of bulgeless disc and normal disc galaxies 
using both the Sersic index and Kormendy relation criteria, 
as described for the main sample. 

Overall, we obtained Sersic parameters in rest-frame B- 
band for 597 bright disc galaxies separated into 211 bulgeless 
disc and 386 normal disc galaxies over three redshift ranges 
(0.77^ 2 <1.0, 0.4^ 2 <0.77 and 0.025C 2 <0.05). For details 
and the catalog please consult fSachdev^ (l2013l l. 


2.2 Obtaining and cleaning the images 

To probe the formation and evolution of bulges in disc galax¬ 
ies with time, it is required to do image analysis and compute 
parameters like concentration and asymmetry of stellar light 
in each individual galaxy of this sample. 

For the main sample (0.45C 2 <1.0), a 10 arcsec cutout 
is downloaded from the HST-AGS data archive and for the 
local sample (0.02^ 2 <0.05), a 3 arcmin cutout is down¬ 
loaded from the NASA-SLOAN ATLAS data archive. Since 
the aim is to do the study in rest-frame B-band, the Hlter 
chosen for obtaining the galaxy image is according to the 
redshift of the galaxy. The images are thus taken in V, i, z 
and g Hlters for redshifts 0.4-0. 6 , 0.6-0. 8 , 0.8-1.0 and 0.02- 
0.05 respectively. The cutout is such that the centre of the 
galaxy (brightest pixel) is at the centre of the image and an 
average-size galaxy covers not more than 60% of the total 
area. 

Out of 597 downloaded images, 27 source images are 
not taken up for analysis. This is because 12 of these galax¬ 
ies are only partially imaged and the rest of the 15 galaxies 
have a highly fragmented light distribution which in part 
appears to be due to multiple overlapping sources. For the 
latter galaxies, it is difficult to determine and study their iso¬ 
lated light distribution. Also, since there is a lack of proper 
structure and there is a degeneracy of bright patches, there 
is no clear area from where the initial central pixel value can 
be selected which is crucial for the computation of radius, 
concentration and asymmetry. Out of these 15 sources, 3 
are in the high redshift range (0.77-1.0) of the bulgeless disc 
category constituting ~1% of this sample. Five are in the 
middle redshift range (0.4-0.77) of the normal disc category 
constituting less than ~4% of this sample. The rest of the 
15 are in the high redshift range (0.77-1.0) of the normal 
disc category and they also constitute ~4% of this sample. 
All these galaxies cover a large range of luminosity, Sersic 
index and half-light radius values. Their small fractions and 
even distribution in terms of redshift and parameter values 
indicate that their removal from the total sample should not 
affect the statistical estimations. 

The major task involved now is to clean or decontami¬ 
nate the 570 galaxy images, i.e., to remove the neighbouring 
sources. These images are final in the sense of flat-fielding, 
bias subtraction, cosmic ray removal etc. To clean the im¬ 
ages, each galaxy image is taken up separately and the value 
of the pixels covered by a neighbouring source is replaced 
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with the average value of the sky pixels surrounding that 
source. The neighbouring sources are recognized through the 
use of SExtractor’s catalog and segmentation map. Along 
with all the individual object coordinates, the SExtractor 
also provides an estimation of the radius of the object that 
contains more than 90% of its light. The replacing of pixels 
is done using IRAF (Image Reduction and Analysis Facil¬ 
ity) imedit task, which creates a circular annulus of a chosen 
radius around the selected central coordinates of the source 
and replaces the pixels within the inner circle with the av¬ 
erage value of the pixels inside the annulus. We have tried 
to ensure that the masking process in each of the galaxy 
images does not strongly affect the light distribution of the 
outer parts of the main source. 

Next important thing is to remove the background flux 
from all the galaxy images. To estimate the background flux, 
a blank patch (i.e. which even without masking is devoid of 
any light sources) is selected near the galaxy. The patch is 
chosen in such a way that it is small and far enough from the 
source, so as not to have any diffuse light, and yet large and 
near enough to give a reasonable estimate. The HST object 
mosaics have an area of 333x333 pixels (10 arcsec squared, 
0.03 arcsec/pixel) and the selected background patch in each 
mosaic covers an area of 30x30 pixels. The SDSS object 
mosaics have an area of 440x440 pixels (3 arcmin squared, 
0.396 arcsec/pixel) and the selected background patch in 
each mosaic covers an area of 40x40 pixels. The mean flux 
per pixel is then estimated from this patch and subtracted 
out from the image. The stability of the process is ascer¬ 
tained by selecting a large number of blank patches for some 
of the galaxy images. It is seen that they provide consistent 
background flux values. After cleaning and background flux 
subtraction, the images are ready for the computation of 
various parameters. 


2.3 Computing Petrosian radii 


The total radius of a disc galaxy is difficult to be deter¬ 
mined in a reproducible manner due to difficulties regarding 
its extent in terms of the dark matter halo. Also, to mea¬ 
sure the optical or stellar extent, use of a hxed isophotal 
value of brightness is not optimal, due to the vast range of 
surface brightnesses in which galaxies exist. However, it is 
important to be determined in such a way that sizes and 
the associated parameters can be compared for galaxies of 

different luminosities and distances. _ 

In that regard, the Petrosian radius (IPetrosianI Il976l l 
is considered an effective way of measuring galaxy sizes 


llBershadv. Jangren fc Conselicell2000l : iGraham et al. l2005ll . 

This radius is determined by tracking the ratio of surface 
brightness at each successively increasing radius to that av¬ 
eraged inside the radius. When the ratio (r;=I(r)/<I(<r)>) 
falls to a chosen small fraction (we take 0.2), that radius 
is multiplied by a factor (we take 1.5) to obtain the Pet¬ 
rosian radius. The values 0.2 and 1.5 are determined to be 
most appropriate (IConselicell2003l : iLotz. Primack fc Madar] 
[20 o 3). For our sample, we first determine the centre accu¬ 
rately using the IRAF task imcentroid and then measure 
the surface brightness at and inside successively increasing 


radii. Then, using 


Rp = l.^*{rp), (15) 

we measure Rp (the Petrosian radius) in rest-frame R-band 
for all the galaxies. This radius is in pixels and converted 
into kpc according to the redshift and chosen cosmology. 


2.4 Computing Concentration 

For the well resolved low redshift galaxies, the bulge promi¬ 
nence is quantified and stud i ed th r ough e.g. bulge/dis c 
decompositions (e.g. iGadottil l2008l : ISimard et al.l li^Olll l. 
However, for high redshifts, single/global Sersic indices 
are commonly used to obtain a quantitative assess¬ 
ment of the bulge component (e.g. iBuitraeo et al.l l2013l : 
iMosleh. Williams fc FraM l2013l . and references therein). 
We have global Sersic index values for all the galaxies in 
our sample, obt ained earlier from Sersic function fitting 
llSachdeval l2013l ~) . We now compute concentration, which 
is considered a more robust pa rameter in terms of surface 
brightness dimming (detailed in iGraham et al.l 2005 ) and a 
better estimator of the bulge-to-total ratio I Gadotti 1200911 . 

To compute concentration, we first measure the to¬ 
tal flux (or total counts) inside the Petrosian radius, i.e. 
the full flux from the source. Then, the number of counts 
in successively increasing radii, from the centre, are com¬ 
puted for each galaxy. When the number of counts is 20% 
of the total number of counts, that radius is taken as 
the ‘inner-radius’. When it’s 80%, that radius is taken as 
the ‘outer-radius’. Concentration of the s ource is defined 
as I Bershndw JaMren fc Conselicel 120001 : IConselit^ l2003l : 
IGraham et al.l 2005l i : 


G = 5.1ogio(^). (16) 

^20% 

Concentration of stellar light, in rest-frame R-band, is 
thus computed for all the galaxies in the sample. The higher 
this measure, larger is the fraction of total light contained 
in the central region. 


2.5 Computing Asymmetry 

Asymmetry in the stellar light of disc galaxies arises from 
features like bars, star forming clumps, spiral arms, rings etc. 
It is also observed to be higher for galaxies which are going 
through interactions or mergers with companion galaxies or 
accreting satellites (or non- virialized baryonic matter) from 
the inter galactic med ium (|Conselicel 1200.31 : iReichard et al.l 
I 2 OO 8 I : llTotz et al.]l2008ll . Since the presence of these features, 
as well as interactions and mergers, are the expected cause 
of bulge formation and evolution, tracking the evolution of 
the asymmetry measure and its relationship with other pa¬ 
rameters is of utmost importance. 

To compute asymmetry of stellar light in 
a disc galaxy, we follow th e pr ocedure given in 
IConselice. Bershadv fc .langreni (l200fll l. We take the 
cutout galaxy image and rotate it around its centre by 180 
degrees. The extraction radius for rotation is given as the 
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Petrosian radius. The rotated image is then subtracted 
from the main image to get the residual image. The flux 
from this residual image is a measure of the flux from the 
asymmetric features of the galaxy. This flux is normalized 
with respect to the total flux from the main image to get 
the asymmetry parameter. 

For the asymmetry parameter to be meaningfully com¬ 
parable for the range of disc gala xies, all the sourc es of prob¬ 
able biases need to be removed l|Conselicel l200.ll 'l . The first 
concern is that the extraction radius should be bias-free, and 
thus, we choose the Petrosian radius. The second concern is 
the centre for rotation, which can produce spurious results 
if not chosen properly. To minimize this effect, the centre 
for rotation is found in an iterative manner: it is the posi¬ 
tion for which the asymmetry of the source attains a global 
minimum. 

Another concern is noise, in terms of the background 
asymmetry. To take that into account, the asymmetry must 
also be computed from an empty background patch and sub¬ 
tracted from the source’s asymmetry. The difficulty with es¬ 
timating the background asymmetry value is that even if 
there is a small amount of diffuse light which is more con¬ 
centrated on one side of the background patch, the asym¬ 
metry value increases by a considerable amount. To avoid 
this issue, we identify relatively clear background patches 
from a number of galaxy images which are chosen such that 
they cover a wide range of RA and Dec. The asymmetry is 
computed on all such patches according to the procedure de¬ 
scribed above. The mean of these values is then used as the 
background asymmetry. This procedure is done separately 
for the galaxy samples in the three redshift ranges. 

The asymmetry of the stellar light distribution of the 
galaxies, in rest-frame R-band, is thus computed for our 
sample. 


2.6 Computing effective stellar mass densities 


The stellar mass density of a galaxy is seen to be cor¬ 
related with its integrated colour, Sersic index, concen¬ 
tration, e nvironmental density, as well as star forma¬ 
tion rate (|KauflmamL_eLjrk 200^ ^;inchmann_et_al. 200j; 


iBaldrv et al.ll2006l : [Driver et al.ll2006l : lBamford et al.l 2009ll . 
The stellar mass density inside the effective radius (termed 
here as effective stellar mass density, ESMD), being highly 
correlated with the bulge properties of the disc galaxy, is a 
useful parameter to be examined. 

To compute this parameter, the first thing is to ob¬ 
tain the stellar mass. For that, we need to multiply total 
luminosity in a given band with the corresponding stellar 
mass-to-light ratio of the galaxy. Total luminosity (L) is ob¬ 
tained in units of solar luminosity from the earlier computed 
absolute magnitudes of the galaxies: 


Mb = -2.5 logio(—)-b 5.38, 

Lq 


(17) 


where Lq is the solar luminosity and 5.38 is solar absolute 
magnitude in rest-frame B-band. The luminosity inside the 
effective radius (or half-light radius) is, by definition, half of 
the total luminosity. 

The next step is to obtain the stellar mass-to-light 
ratio which is known to strongly co rrelate with the inte¬ 
grated colours jBell fc de Jond 200^. To obtain the rest- 


frame colours, we employ EAZY, which is a photometric red- 
shift c ode that provides estimates for rest-fra me colour in¬ 
dices (IBrammer, van Dokkum fc CoPDill2008ll . It compares 
photometric data to the synthetic photometry of a large 
range of template spectra and outputs the best match. The 
important feature of this code is that if the redshift of 
the source is known with reasonable accuracy, it can be 
given as a prior or held fixed. Also, the template spectra 
are based on semi-analytical models and not on empiri- 
cal spectroscopic samples, which are us ually highly biased 
dBrammer. van Dokkum fc CoppI l2008l l. Note that there 
are o ther codes similar to EAZY, e.g., FAST dKriek et al.l 
I2OO9II and LEPHARE dArnouts fc Ilberf^l2011^ . 


The stellar mass-to-light ratio is then computed 
from the rest- f rame B — V colour, using values from 
iBell fc de Jond d200lll for the mass-depend ent galaxy for¬ 
mation model with bursts dCole et al.ll2000f) . We chose this 
model because mass dependence is the common fe ature of 
the p resently acceptable galaxy form ation scenarios I Benson] 
l2010l . and references therein). Also. [Bell fc de Jong 1 2001 ) 
selected this model as their default model claiming that this 
model reproduces the trends in age and metallicity with re¬ 
spect to surface brightness with the least scatter for local 
spiral galaxies. 

We, thus, obtain the total stellar mass for each galaxy 
in our sample using the stellar mass-to-light ratio computed 
above. This method of obtaining stellar masses is extensively 
used in extra-galactic studies where well resolved spectral 
data is not available. The concern is that we are not taking 
the galaxy colour gradient into account by opting for its 
global colour. This may lead to an underestimation of mass, 
however, it should not affect our analysis of relative increase 
in mass and density. 

The stellar mass inside effective radius is half of the 
total stellar mass. This mass is then divided by the area 
within the effective radius to obtain stellar mass density 
inside the effective radius (or ESMD) in units of solar mass 
per square kpc. 


2.7 Checking for accuracy and error computation 

The working of the overall procedure/code written to com¬ 
pute Petrosian radii, concentrations and asymmetries of the 
sources was tested using artificially created images and real 
images with known parameters. 

The artificial images are created using Galflt’s Sersic 
component. First, we fix the Sersic index (at n=l) and cre¬ 
ate images with varying half-light radius, i.e. re=40, 50, 60, 
62, 65 and 70 pixels. Since the Sersic index is fixed, the 
Petrosian radius is expected to increase with the increase 
in half-light radius. This is indeed seen as the Petrosian 
radius for these images is computed to be 98.46(±4.92), 
122.74(±6.14), 146.65(±7.33), 151.47(±7.57), 158.44(±7.92) 
and 169.91 (±8.49) pixels respectively. 

Similarly, for fixed half-light radius (at re=50 pixels), 
we create images with varying Sersic index, i.e. n=0.8, 1.2, 
1.6, 2.0 and 2.4. Here, since the radius containing half of the 
light is held fixed, the concentration of the galaxy is expected 
to increase with the increase in its Sersic index. This is in¬ 
deed reported as the concentration index for these images 
is computed to be 2.65(±0.07), 3.01(±0.08), 3.29(±0.09), 
3.52(±0.09) and 3.71(±0.10) respectively. 
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Next, we create images with different levels of asymme¬ 
try. This is done by keeping all the parameters (apparent- 
magnitude, half-light-radius, Sersic index) fixed and adding 
Fourier modes. A detailed discussion of the representation 
of various asy mmetric gal a ctic fe atures with Fourier modes 
is provided in IPeng et al.l (|201(J) . We add the first Fourier 
mode of varying amplitudes, i.e. 0.07, 0.09, 0.11, 0.13, 0.15, 
0.17 and 0.19. The asymmetry measure responds favourably 
such that it increases with the increasing amount of the 
Fourier mode amplitude. It is computed to be 0.080, 0.103, 
0.125, 0.147, 0.169, 0.191 and 0.213 respectively with error 
in the range of ±0.045 to 0.048. The procedure, thus, repro¬ 
duces the measures with reasonable accuracy, responding 
well to the small shifts produced in the structural proper¬ 
ties. 

A concern relating to the accurate computation of 
these parameters at high redshift is the cosmological surface 
brightness dimming, which may lead to the non detectabil¬ 
ity of faint features. However, the surface bri ghtness evolu¬ 
tion, reported to be 1-1.5 mag since z = 1 ll Barden et al.l 
I 2 OO 5 I : [Melbourne et al.ll2007l : ISachdev^l2013l i is expected to 
counter this dimming. Another concern is that of low reso¬ 
lution of galaxies at high redshift. It should not affect our 
study (till 2 ~ 1) because the images are from overlap¬ 
ping 5-orbit depth GOODS survey using HST-ACS, which 
provides a resolution of 0.03 arcsec per pixel . This is veri¬ 
fied in a quantitative way by IConselicd (l2003l i . They simu¬ 
lated the local bright galaxies to higher redshifts as per how 
these galaxies would be imaged by various surveys. Then 
they compared the values measured at « 0 to the val¬ 

ues measured at various redshifts. For z = 1 (for GOODS 
HST ACS), they report a marginal change of 0.10T0.18 in 
the concentration measure and a change of -0.03T0.07 in 
the asymmetry measure. Thus, the indices are highly repro¬ 
ducible with negligible scatter. 

The error bars in the measurement of these parame¬ 
ters stem from the uncertainties involved in determining the 
total flux associated with the pixels of interest and also in 
the selection of the pixels of interest. The error associated 
with the measurement of flux using the aperture photometry 
package apphot of IRAF is calculated using: 

error = v(counts/gain ± area * stdev^ 

2 2 (18) 

± area * stdev /nsky), 

where gain is in electrons per adu and area (of the aperture) 
in square pixels, stdev is the standard deviation of the sky 
counts and nsky is the number of sky pixels. 

Out of the total sample of 570 galaxies, the algo¬ 
rithm/procedure did not converge to give the parameter val¬ 
ues for 3 galaxies, reporting a floating point zero error. This 
is an extreme case caused for sources with central point of 
such high brightness that the inner radius goes to zero. 


2.8 Overall data sample 

We obtained parameters for overall 567 bright (Ms ^-20) 
bulgeless (i.e. without classical bulge) and normal (i.e. with 
classical bulge) disc dominated galaxies in the three red¬ 
shift ranges (263 in 0.77^ 2 : <1.0, 203 in 0.4^ z <0.77, 
101 in 0.02^ z <0.05) in rest-frame R-band. We have their 
redshifts, absolute magnitudes, half-light radii in kpc and 
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Figure 1. Some of the bulgeless disc galaxies lying in the three 
redshift ranges (0.02-0.05, 0.4-0.77 and 0.77-1.0) are shown. The 
non-parametric measures computed for each source are shown at 
the top left corner. They are in this particular order: ID (as per the 
NYU-VAGC and GOODS HST-ACS catalogs), Petrosian radius. 
Concentration and Asymmetry. Precise positional information of 
the galaxies is provided in Table [T] HST-ACS galaxies cover ~7 
arcsec (out of the 10 arcsec image cutout) and SDSS galaxy im¬ 
ages cover ~2 arcmin (out of the 3 arcmin image cutout). 


Sersic index from ISachdeval (l2013l ). We have now computed 
their Petrosian radii in kpc, concentration, asymmetry, total 
stellar mass and effective stellar mass density. Some of the 
bulgeless and normal disc galaxy images from the sample 
are shown in Fig. [T] and Fig. [2] for the three redshift ranges. 
The computed parameters along with the associated errors 
for these images are provided in Table [T] In the next section, 
we present the results that provide insights into bulge forma¬ 
tion and evolution occurring in disc galaxies, by examining 
the evolution and relationships of these parameters. 
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Table 1. Parameters for galaxies shown in Fig. [J and Fig. [2] 


Source 

ID 

RA (J2000) 
degrees 

Dec (J2000) 
degrees 

Z 

C 

C — err 

A 

A — err 

Rp 

Kpc 

Rp — err 
Kpc 

SDSS 

107284 

18.182632 

15.70782 

0.039 

2.318 

0.118 

0.492 

0.026 

20.962 

1.048 

SDSS 

138496 

140.67966 

57.51786 

0.049 

2.106 

0.128 

0.359 

0.036 

21.709 

1.085 

SDSS 

236819 

318.44444 

-5.8170589 

0.049 

2.543 

0.157 

0.367 

0.034 

21.686 

1.084 

SDSS 

364153 

148.98351 

53.693652 

0.044 

2.531 

0.127 

0.356 

0.027 

22.805 

1.140 

SDSS 

63945 

226.72858 

0.1863334 

0.035 

3.586 

0.182 

0.222 

0.016 

18.186 

0.909 

SDSS 

87659 

231.89702 

0.28918961 

0.049 

3.333 

0.254 

0.192 

0.009 

16.739 

0.837 

SDSS 

107614 

28.966593 

14.940272 

0.044 

3.042 

0.203 

0.397 

0.018 

17.078 

0.854 

SDSS 

110067 

42.326159 

-8.1749502 

0.030 

3.485 

0.156 

0.211 

0.008 

18.788 

0.939 

HST-ACS 

1355 

53.0167785 

-27.7189942 

0.538 

2.93 

0.192 

0.376 

0.034 

9.661 

0.483 

HST-ACS 

14303 

53.1091768 

-27.8529078 

0.629 

2.937 

0.096 

0.759 

0.119 

21.586 

1.079 

HST-ACS 

28155 

53.190166 

-27.7349438 

0.530 

2.417 

0.104 

0.643 

0.077 

15.897 

0.795 

HST-ACS 

31666 

53.2226372 

-27.8475767 

0.432 

2.755 

0.121 

0.598 

0.058 

12.901 

0.645 

HST-ACS 

2587 

53.0311719 

-27.7357189 

0.935 

2.483 

0.209 

0.937 

0.197 

10.912 

0.546 

HST-ACS 

11609 

53.093944 

-27.8727229 

0.789 

2.573 

0.202 

0.581 

0.074 

9.875 

0.494 

HST-ACS 

19859 

53.1384525 

-27.6806527 

0.828 

2.581 

0.207 

0.657 

0.111 

9.887 

0.494 

HST-ACS 

27617 

53.1859173 

-27.7756062 

0.848 

2.155 

0.186 

0.791 

0.16 

9.571 

0.478 

HST-ACS 

3467 

53.0391609 

-27.7100294 

0.716 

3.12 

0.276 

0.312 

0.036 

8.484 

0.424 

HST-ACS 

12196 

53.0975165 

-27.7212663 

0.551 

3.09 

0.156 

0.531 

0.048 

13.391 

0.669 

HST-ACS 

19467 

53.1361845 

-27.836454 

0.633 

3.466 

0.139 

0.46 

0.05 

17.084 

0.854 

HST-ACS 

26780 

53.179339 

-27.9235978 

0.615 

3.296 

0.151 

0.596 

0.07 

15.290 

0.764 

HST-ACS 

4982 

53.050917 

-27.7724075 

0.924 

3.211 

0.246 

0.482 

0.043 

10.156 

0.508 

HST-ACS 

9837 

53.0831774 

-27.7471694 

0.883 

2.937 

0.163 

0.678 

0.108 

14.526 

0.726 

HST-ACS 

10198 

53.0854551 

-27.6830331 

0.847 

3.224 

0.138 

0.837 

0.129 

17.162 

0.858 

HST-ACS 

16603 

53.120832 

-27.8230569 

0.783 

3.067 

0.367 

0.272 

0.033 

6.149 

0.307 


3 RESULTS 

The formation and evolution of bulges in disc galaxies can 
be probed by examining the mutual evolution and relation¬ 
ship of those parameters whose measure is associated with 
bulge properties. We, thus, examine the evolution of size, 
concentration, inner stellar density and asymmetry. 


3.1 Size evolution 


The sizes of the massive disc galaxies in the universe 
are seen to undergo a dramatic increase with time (e.g. 



UuhllQ_et_alj2QQ^ Buitra go et al.ll2008l: Ivan Dokkum et 


Carrasco. Conselice fc Trujillo! 20101 : ICassata et al 


Some studies have attempted to explain this in¬ 
crease through inside-out processes, m ajor mergers, minor 
mergers, accretion. AGN processes etc. (iHookins et ]l200fll: 
Kavh;aj_et_^lJ_ 2001 ; Bluck et al.l I 2 OI 2 I : lOwnsworth et all 


2012l:lMcLure et al.ll201^ 


The mechanisms that lead to the overall increase in 
disc length are also expected to produce variations in the 
inner region. For example, secular evolution driven by var¬ 
ious asymmetric structures causes galactic discs to expand 
on the outside and cont r act on the inside llTremainelll989l : 
Ivan Dokkum et al ] l20inl: [k ormen d^ l201.‘ll '). To understand 
the relative role of the mechanisms in disc growth, the in¬ 
crease in the total optical extent should be tracked with the 
change in the inner region properties. 

We, therefore, examine the increase in the Petrosian 
radius of the galaxy relative to the increase in its half-light 
radius. This relative increase will be seen separately for discs 
with and without classical bulges at the three redshift ranges 
(0.02-0.05, 0.4-0.77, 0.77-1.0). 


The distribution of the Petrosian radius for bulgeless 
and normal disc galaxies is shown in Fig. [3] for the three 
redshift ranges. For the bulgeless disc galaxies it increases 
by 85(±4)% (from 12.71 kpc to 23.46 kpc) from 2 - 0.9 
to the present epoch. Over the same time range, for the 
normal disc sample, it increases by 94(±6)% (from 10.50 
kpc to 20.37 kpc) (see Table [2J. 

Thus, both morphological types show a significant in¬ 
crease in size. However, the evolution is faster and leads to 
a more effective increase in size for the normal disc galax¬ 
ies. Nevertheless, it is interesting to note that bulgeless disc 
galaxies are larger than the normal disc galaxies, on average, 
at all redshift bins. 

In contrast to the Petrosian radius, the half-light radius 
for both bulgeless and normal disc samples shows a minor in¬ 
crease of 28(±2)% (from 5.43 kpc to 6.97 kpc) and 22(±1)% 
(from 3.72 kpc to 4.55 kpc), respectively, from 2 ~ 0.9 to 
the present epoch (Table [21). 

If we see the relative quantity, i.e. the ratio of Petrosian 
radius to half-light radius, it increases from 2.34 to 3.36 from 
2 ~ 0.9 to the present epoch for bulgeless disc galaxies. 
However, the increase is larger, from 2.82 to 4.48, for the 
normal disc galaxies. For these galaxies, their total radius 
becomes ~4.5 times their half-light radius, as we reach the 
present epoch. 

The two sizes (half-light radius and Petrosian radius), 
though computed in totally different ways (parametrically 
and non-parametrically), are found to be highly correlated 
in our galaxy sample, as seen in Fig. [4] The linear fit relation 
found for the highest redshift range (0.77-1.0) almost over¬ 
laps the relation found for the intermediate redshift range 
(0.4-0.77). However, the slope changes drastically for the lo- 
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Figure 2. Some of the normal disc galaxies lying in the three 
redshift ranges (0.02-0.05, 0.4-0.77 and 0.77-1.0) are shown. The 
non-parametric measures computed for the source are shown at 
the top left corner. They are in this particular order: ID (as per the 
NYU-VAGC and GOODS HST-AGS catalogs), Petrosian radius. 
Concentration and Asymmetry. Precise positional information of 
the galaxies is provided in Table[T] HST-AGS galaxy images cover 
~7 arcsec (out of the 10 arcsec image cutout) and SDSS galaxy 
images cover ~2 arcmin (out of the 3 arcmin image cutout). 


cal redshift range (0.02-0.05): 

Rp = 1.69(±0.06) * + 3.62(±0.31); 0.4 < z < 1.0, (19) 

Rp = 4.92(±0.68) * - 2.04(±0.59); 0.02 ^ 2 : < 0.05. (20) 

The change of slope further signifies that the half-light ra¬ 
dius reduces to a much smaller fraction of the Petrosian 
radius as we reach the present epoch. 

Overall, we find that while the total extent of our discs 
almost doubles with time, the radius containing half the stel¬ 
lar light increases marginally in comparison. This suggests 


significant peripheral increase, which is seen to be somewhat 
more pronounced for the normal disc sample. 


3.2 Concentration evolution 


The Sersic index is extensively used not just for morpho¬ 
logical characterizati on of galaxies but also to study their 
structural evolut ion dShen et al.ll2003l: iBlanton et ahl 2005^ 


Fislyr^jDrorv _ 20od;_ Conselice et al. 2011 : Bruce et al.l 


2 OI 2 I : Buitrago et al.ll2013l ). Being related to the steepness of 


the intensity profil e, it is a measure of the promin ence of the 
bulge component llSersid[l968l : IPeng et al.ll2002l) . The con¬ 
centration parameter, however, is seen to be a bett er estima¬ 
tor o f the bulge presence for low redshift galaxies (IConselicd 
I 2 OO 3 I : [Graham et al.ll2005l : lGadotBll2009l ). 


Our full (disc dominated) sample is separated on the 
basis of Sersic index (along with the Kormendy relation). 
We now examine the evolution in their concentration value. 
The distribution of the concentration for bulgeless and nor¬ 
mal disc galaxies is studied for the three redshift ranges 
(0.02-0.05, 0.4-0.77, 0.77-1.0). The mean concentration of 
the bulgeless disc galaxies shows a statistically insignificant 
increase from 2 ~ 0.9 to the present epoch (Table EJ. How¬ 
ever, for the normal disc sample it increases by 12.3(±0.3)% 
(from 2.77 to 3.11), over the same time range (Table[3]). 

Thus, while the average bulgeless disc galaxy concentra¬ 
tion remains almost similar over the three redshift ranges, 
there is a significant (albeit admittedly small) increase in 
average normal disc galaxy concentration. 

While the Sersic index has been obtained parametri¬ 
cally, i.e., by fitting a function to the surface brightness pro¬ 
file, the concentration is based on total count ratios, i.e. 
non-parametrically. We examine the relationship of the two 
parameters for the full sample (Fig. [5]) in the three redshift 
ranges. The two quantities are seen to be well correlated 
for both morphological types and follow a single relation for 
both the higher (0.77-1.0) and intermediate (0.4-0.77) red¬ 
shift ranges: 


C = 0.38(±0.03) *n + 2.32(±0.04), 0.4 < 2 < 1.0 (21) 


The two quantities appear to provide a similar estimate of 
the bulge component. Thus, according to our study they 
are equally well deserving to be chosen at high redshifts for 
morphological determination. 

However, the correlation is absent for the local redshift- 
range (0.02-0.05). This lack of correlation at local redshifts 
is also reported by Gadotti (2009). We speculate that for 
highly resolved local galaxies the intensity gradient between 
the bulge and the disc is enhanced, leading to poorer fits 
when only a single function is used to fit the entire galaxy. 


3.3 Stellar density evolution 


Total stellar mass is one of the most significant properties 
of a galaxy and is seen to be correlated with not just the 
overall conc entration but also the star formation rate of 
the galaxy llGaon. Capaccioli fc D’Onof^ Il993l: iGonselice 


2002; iNoeske et al. 2007 : Disney et al. 20081 : iBauer et al 
201111 . Recent studies have found that close to half of 
the p resent stellar mass of the galaxies assembled by 
2 ~ 1 llBundv. Ellis fc Gonselic^l200d : [Mortlock et al.ll201li : 
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Figure 3. Distributions of Petrosian radius for bulgeless disc (solid lines) and normal disc (dashed lines) galaxies for the three redshift 
ranges. The distribution of the means of the two samples with redshift is also shown. There is a significant increase in sizes for both 
samples with time. 




Figure 4. Petrosian radius is plotted against half-light radius for bulgeless and normal disc galaxies in the three redshift ranges. The 
solid line in each graph marks the linear relation followed by the full sample at that redshift range. The highest redshift range (0.77-1.0) 
relation is shown on the two lower redshift range (0.4-0.77, 0.02-0.05) plots with dotted lines. It is seen to almost overlap the intermediate 
redshift range (0.4-0.77) relation. However, for the local sample (0.02-0.05), the slope is entirely different. There the half-light radius is 
seen to be a much smaller fraction of the Petrosian radius. For the Petrosian-radius, the error scales with the value and is seen to be in 
the range of ±5%. For half light radius, the typical error is of ±0.2 kpc. 
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Table 2. Mean and median values of Petrosian radius and half-light radius. 


Redshift 

range 

Disc 

type 

No. of 

sources 

Mean of 

Petrosian-radius 
< Rp >(Kpc) 

Std.dev. 
of Rp 

cr 

Median 
of Rp 

Mean of 
half-light radius 
< Re >(Kpc) 

Std.dev. 
of Re 

a 

Median 
of Re 

0 . 77 - 1.0 

bulgeless 

105 

12.706(±0.317) 

3.246 

12.249(±0.397) 

5.429(±0.179) 

1.844 

4.902(±0.224) 

0 . 77 - 1.0 

normal 

158 

10.498(±0.325) 

4.087 

9.411(±0.407) 

3.725(±0.149) 

1.885 

3.170(±0.187) 

0 . 4 - 0.77 

bulgeless 

73 

13.787(±0.459) 

3.922 

13.871(±0.575) 

5.947(±0.251) 

2.148 

5.539(±0.314) 

0 . 4 - 0.77 

normal 

130 

11.075(±0.368) 

4.194 

10.507(±0.461) 

4.446(±0.193) 

2.208 

4.080(±0.242) 

0.02-0.05 

bulgeless 

25 

23.464(±1.035) 

5.177 

22.395(±1.297) 

6.972(±0.500) 

2.504 

6.508(±0.626) 

0.02-0.05 

normal 

76 

20.373(±1.119) 

9.761 

18.244(±1.402) 

4.552(±0.147) 

1.280 

4.368(±0.184) 



Figure 5. Concentration is plotted against Sersic index, n, for bulgeless and normal disc galaxies in the three redshift ranges. The solid 
line in each graph marks the linear relation followed by the full sample at that redshift range. The highest redshift range (0.77-1.0) 
relation is shown on the two lower redshift range (0.4-0.77, 0.02-0.05) plots with dotted lines. The two parameters are seen to be highly 
correlated and the relation almost overlaps for the two high redshift ranges (0.77-1.0, 0.4-0.77). However, the local sample shows a lack 
of correlation in the two values. For concentration, typical error on the value is ±0.17 and for Sersic index it is ±0.05. 


iMarchesini et al1l2014l : [Ownsworth et al.]|2014h . We first ex¬ 
amine the growth in the stellar mass for our full sample from 
z ~ 0.9 to the present epoch (Fig. |6]). 

The mean values of the log of total stellar mass (in units 
of solar mass, Mq) are given in the three redshift ranges in 
Table 21 The bulgeless disc galaxies witness an increase of 
3.1x10'^° Mq (from 2.7 to 5.8x10^° Mq) from 2 ; ~ 0.9 to 
the present epoch. Over the same time range, the normal 
disc galaxies witness an increase of 4.5x10^° Mq (from 4.4 
to 8.9x10^° Mq). 

Thus, we hnd that both the morphological types have 
gained more than 50(±6)% of their present stellar mass since 
2 ~ 0.9. The interesting part is that the increase in the total 
stellar mass of a normal disc galaxy is ~1.5 times more than 
that seen for a bulgeless disc galaxy over the last ~8 Gyrs, 
on average. The difference in the average total stellar mass 
of the two morphological types almost doubles (from 1.7 to 
S.lxlO*^® Mq) from 2 ~ 0.9 to the present epoch. 

Next we examine the growth of stellar mass density in 
the inner region i.e. the effective stellar mass density. The 
distribution of the effective stellar mass density on the log- 
scale is shown in Fig.Qfor bulgeless and normal disc galaxies 
in the three redshift ranges. 

Examining the mean values from 2 ~ 0.9 to the present 
epoch (Table HI, we find an increase of 4.7x10^ Mg/kpc^ 
(from 1.6 to 2.1x10® MQ/kpc^), on average, for the effective 


stellar mass density of bulgeless disc galaxies; and of 8.8 x 10^ 
Mg/kpc^ (from 6.6 to 7.4x10® Mg/kpc^), on average, for 
the normal disc galaxies. 

While the normal disc galaxies witness an increase of 
13(±5)% in their effective stellar mass density, this is more 
prominent for the bulgeless disc galaxies, namely 30(±1)%. 
However, in absolute terms, the increase for the average nor¬ 
mal disc galaxy is ~1.8 times that seen for a bulgeless disc 
galaxy over the last Gyrs. Thus, for both total stellar 
mass and effective stellar mass density, the increase seen for 
the normal disc sample is considerably more than that seen 
for the bulgeless disc sample. 

In addition, we note that the effective stellar mass den¬ 
sity can be a better indicator of galaxy morphology than 
concentration. This is seen using the Kolmogorov-Smirnov 
test which is considered an efficient mathematical tool for 
determining the significance of the difference of two distri¬ 
butions. We use this test to compare the bulgeless and nor¬ 
mal disc galaxy samples with respect to concentration and 
effective-mass-density in three redshift ranges (Table O. 

The null hypothesis is that the two distributions are 
from the same parent set and this test quantifies a proba¬ 
bility for this null hypothesis. For that, we compute the dis¬ 
tance (D-observed) between the empirical cumulative distri¬ 
bution functions of the two samples. Based on this distance 
and the sample size, the probability for the null hypothesis 
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Table 3. Mean and median values of concentration and asymmetry. 


Redshift 

range 

Disc 

type 

No. of 

sources 

Mean of 

concentration 

<0 

Std.dev. 
of C 

(T 

Median 
of C 

Mean of 
asymmetry 
< A > 

Std.dev. 
of A 

a 

Median 

of A 

0 . 77 - 1.0 

bulgeless 

105 

2.748(±0.027) 

0.278 

2.772(±0.034) 

0.827(40.014) 

0.145 

0.831(40.017) 

0 . 77 - 1.0 

normal 

158 

2.767(±0.030) 

0.382 

2.784(±0.037) 

0.568(40.014) 

0.181 

0.559(40.017) 

0 . 4 - 0.77 

bulgeless 

73 

2.726(±0.037) 

0.317 

2.677(±0.046) 

0.719(40.019) 

0.159 

0.718(40.024) 

0 . 4 - 0.77 

normal 

130 

2.819(±0.034) 

0.391 

2.776(±0.042) 

0.515(40.016) 

0.188 

0.494(40.020) 

0.02-0.05 

bulgeless 

25 

2.791(±0.067) 

0.448 

2.674(±0.084) 

0.412(40.023) 

0.115 

0.394(40.029) 

0.02-0.05 

normal 

76 

3.107{±0.064) 

0.560 

3.010(40.080) 

0.308(40.015) 

0.131 

0.280(40.019) 


Table 5. KS-test D-observed and Probability values 


Parameter 

Redshift 

range 

D-observed Probability 

Cone 

0.02-0.05 

0.337 

0.021 

Log-EMD 

0.02-0.05 

0.709 

3.4e-09 

Cone 

0.40-0.77 

0.178 

0.092 

Log-EMD 

0.40-0.77 

0.395 

9.7e-07 

Cone 

0.77-1.00 

0.132 

0.204 

Log-EMD 

0.77-1.00 

0.414 

8.9e-10 


is found (Table [S]). In the case of concentration, this proba¬ 
bility is not convincingly low for the two distribntions to be 
considered significantly different. However, it is found to be 
negligible at all redshift ranges in the case of effective stellar 
mass density. 

3.4 Correlations with asymmetry 


The asymmetric features in the disc galaxy, cre¬ 
ated by both internally and externally driven pro¬ 
cesses, lead to the formation of pseudo- and clas- 


sical bulges in disc galaxies ( 

Khochfar & SillJ 

2006|; 

Elmeereen. Bournaud & Elmesreen 

I 2 OO 81 : 1 Joeee et al.l 

200 g|; 

Kormendv et al.ll2010l: Honkins et al.ll201C 

: Conselicd 20 14). 


Thus, it is imperative to analyse the evolution of asymmetry 
and its relationship with other parameters of disc galaxies. 

The distribution of the asymmetry for bulgeless and 
normal disc galaxies is shown in Fig.|S]for the three redshift 
ranges. The mean asymmetry for bulgeless disc galaxies falls 
by 50(±3)% (from 0.83 to 0.41) from 2 ~ 0.9 to the present 
epoch (Table|3J. Over the same time range, the mean asym¬ 
metry for the normal disc sample falls by 45(±2)% (from 
0.57 to 0.31; Table E]). 

Both bulgeless and normal disc galaxies show a huge de¬ 
cline in their asymmetry value with time, indicating the dis¬ 
appearance of asymmetric features. The bulgeless disc sam¬ 
ple is more asymmetric than the normal disc sample at all 
redshift ranges (Fig.|8}. However, due to the more significant 
fall seen in the average asymmetry value of the bulgeless disc 
sample, it reaches closer to the average asymmetry value of 
the normal disc sample. 

The asymmetry parameter shows a correlation with 
half-light radius (Fig. El). Bulgeless disc galaxies, having 
larger half-light radii than the normal disc galaxies, have 
higher asymmetries. At fixed radii, bulgeless disc galaxies 


are found to be more asymmetric than normal disc galaxies, 
on average. 

The asymmetry parameter shows an anti-correlation 
with effective stellar mass density, as explored in Fig. [10] 
for both morphological types. There is a steep fall in the 
asymmetry with the increase in effective stellar mass den¬ 
sity for all the redshift ranges. The slope for 0.4^ 2 <1.0 is: 

A = -0.21(±0.02) * log(EMD) + 2.47(±0.13). (22) 

Bulgeless disc galaxies being more asymmetric and less 
dense as compared to normal disc galaxies are on the higher 
end of the slope in all three redshift ranges. 


4 DISCUSSION 

We have a total sample of 567 disc-dominated galaxies sep¬ 
arated into bulgeless (without classical bulge) discs and 
normal (with classical bulge) discs in three redshift ranges 
(0.77^ 2 <1.0, 0.4^ 2 <0.77 and 0.02< 2 <0.05). We have 
examined the evolution in size, concentration, effective stel¬ 
lar mass density and asymmetry of these two samples. We 
first list our major findings and then discuss the implica¬ 
tions: 

(i) Both morphological types show a significant increase 
[94(±6)% for the normal disc sample] in total optical extent 
since 2 ~ 0.9. The half-light radius of the galaxies witnesses 
a much smaller increase [ 22 (± 1 )% for the normal disc sam¬ 
ple] in comparison. This peripheral size evolution is more 
evident for the normal disc sample in which the outer ra¬ 
dius becomes ~5 times the inner radius by 2 ~ 0 . 

(ii) The mean concentration of stellar light undergoes a 
significant increase (AC = 0.34) for the normal disc sample 
as compared to the bulgeless disc sample (AC = 0.04) since 
2 ~ 0.9. The concentration parameter of the galaxy is seen to 
be well correlated with its global Sersic index value for both 
morphological types at higher redshift ranges (0.44 — 0.77, 
0.77 — 1.0). However, the correlation is absent for the local 
sample (0.02 — 0.05). 

(iii) Both morphological types have gained more than half 
of their present stellar mass since 2 ~ 0.9. In absolute terms, 
the increase in the total stellar mass as well as effective stel¬ 
lar mass density is significantly more important (~1.5 and 
1.8 times respectively) for the normal disc sample. 

(iv) The bulgeless disc sample is more asymmetric than 
the normal disc sample at all redshift ranges. Both samples 
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Figure 6. Distribution of log of total stellar mass (in units of Mq) is shown in the three redshift ranges for bulgeless disc (solid lines) 
and normal disc (dashed lines) galaxies. The distribution of their mean values with redshift is also shown. Since it is on log scale, even 
a slight shift corresponds to a huge increase in the mass of the galaxy. A shift towards higher mass can be seen for both morphological 
types. Their stellar mass at 2 : ~ 0 is more than double of that at 2 ~ 0.9. 


Table 4. Mean values of log of total stellar mass (LTSM) and log of effective-stellar-mass-density {LESMD). 


Redshift 

range 

Disc 

type 

No. of 

sources 

Mean of 

LTSM 
< LTSM > 

Std.dev. 
of LTSM 

a 

Median 

of LTSM 

Mean of 
LESMD 
< LESMD > 

Std.dev. 
of LESMD 

(7 

Median 

of LESMD 

0 . 77 - 1.0 

bulgeless 

101 

10.438(±0.057) 

0.574 

10.381(±0.071) 

8.214(±0.058) 

0.584 

8.219(±0.073) 

0 . 77 - 1.0 

normal 

151 

10.647(±0.048) 

0.593 

10.610(±0.060) 

8.817(±0.049) 

0.598 

8.795(±0.061) 

0 . 4 - 0.77 

bulgeless 

71 

10.503(±0.072) 

0.603 

10.386(±0.090) 

8.211(±0.064) 

0.544 

8.178(±0.080) 

0 . 4 - 0.77 

normal 

121 

10.634(±0.055) 

0.609 

10.562(±0.069) 

8.654(±0.062) 

0.682 

8.575(±0.078) 

0.02-0.05 

bulgeless 

25 

10.768(±0.056) 

0.280 

10.755(±0.070) 

8.324(±0.049) 

0.248 

8.339(±0.061) 

0.02-0.05 

normal 

76 

10.950(±0.045) 

0.391 

10.943(±0.056) 

8.872(±0.043) 

0.377 

8.854(±0.054) 


witness a fall in their mean asymmetry value from 2 ~ 0.9 to 
2^0, the fall being more drastic [~50(±3)%] for the bulge¬ 
less disc sample. Asymmetry is found to be correlated with 
the half-light radius, and anti-correlated with the effective 
stellar mass density of the galaxy. 


4.1 Impact of internal evolution 

In our sample, from 2 0.9 to 2 ~ 0, the Petrosian ra¬ 

dius increases more significantly than the effective radius. 


In other words, the strong increase in the total radius of 
the full sample is not reflected in its half-light radius. This 
peripheral increase appears to provide evidence in support 
of internal secular evolution in which, due to the outward 
transfer of angular momentum, galaxy disks are expected to 
expa nd on the outside and get more concentrated on the in- 

side (lTremainelll989l:ICombesll200ll:[Kormendv &: KennicuttI 

l2004h . However, there are complexities with this explana¬ 
tion. 


The asymmetric features such as spiral arms and bars, 







































































14 


Sonali Sachdeva, Dimitri A. Gadotti, Kanak Saha, Harinder P. Singh 





9.25 - 


A 

o' 

s 

CO 

W 


8.75 - 


ty) 8.5 
O 

V 

8.25 


« bulgeless 
o normal 


0 0.5 1 

Redshift 


Figure 7. Distribution of log of effective stellar mass density is shown in the three redshift ranges for bulgeless disc (solid lines) and 
normal disc (dashed lines) galaxies. The distribution of their mean values with redshift is also shown. Since it is on log scale, even a 
slight shift corresponds to a huge increase in the mass density of the galaxy. The normal sample shows fluctuation, however, the density 
at 2 : ~ 0 is more than that at 2 : ~ 0.9 for both morphological types. 


that induce and speed up internal secular evolution llJoeee 

2ooe 

: Kormendv & Kennicutt 2004 

; Coelho & Gadotti 

2011 

: ISheth et al.ll2012l: ICheune et al. 

|2013|) are known to 


be p resent in th e bulgeless disc galaxies with greater propen¬ 
sity u t3l2013ll . In our study also the bulgeless disc sample 
is found to be more asymmetric than the normal disc sam¬ 
ple at all three redshift ranges. This in turn might favour 
the build-up of central concentrati on and eventual fading of 
the asymmetric features w ith time llKormendv fc KennicuttI 
l2004 lAth anassoui3l2005l i in bulgeless galaxies. This can be 
the probable cause for a relatively steeper decrease in the 
asymmetry value of our bulgeless disc sample. 

Naively, thus, internal secular evolution is expected to 
be more efficient for the bulgeless disc sample. However, we 
find that the evolution in terms of size and increase in the 
density of the inner-region, as seen through concentration 
and effective stellar mass density, in absolute terms, is con¬ 
siderably more for the normal disc sample. This gives us an 
indication that internal secular processes are not the only 
evolution determining forces. 

In addition to that, it is known through simulations 
and observations that the huge increase in disc galaxies’ 
total optical extent is occurring due to s tellar mass build 
up in the outer regions o f these galaxies lICaPDellari et all 
I2OO9I : iNewman et al.|[2010l l. We also report growth in stellar 


masses by a factor of two since z~l for both the morpho¬ 
logical samples. 

This increase in size in terms of stellar mass build¬ 
up as well as the increase in internal density cannot be 
achieved only through the rearrangement of mass and 
angular-momentum. Thus, in addition to the internal sec¬ 
ular evolution, there might be other processes causing the 
disc evolution from z~0.9 to the present epoch. We argue 
that such inside-out growth might be driven by the trans¬ 
fer of matter from the circumgalactic environment to the 
galaxy, in the next section. 


4.2 Impact of external evolution 


Our findings that discs have grown by such large factors 
from 2 : ~ 1 to now argues that stellar discs are robust struc¬ 
tures, difficult to be destroyed, and that catastrophic merg¬ 
ers are rather rare at the second half of the age of the uni¬ 
verse. This is in agreement with recent studies, based on 
simulations and observations, which repo rt a signihcant de¬ 
cline in the major merger rate with tim e llJoge e et al.ll2009l : 
IConselice. Yang fc Bluckll2009l : iBluck et aLl2012h . 

However, there is continuous accretion of matter 


from the inter galactic medium, and minor mergers are 
also frequent IParrv. Eke fc FrenkI l2009l : iKavira I I2OIOI: 
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Figure 8 . Distribution of asymmetry for bulgeless disc (solid lines) and normal disc (dashed lines) galaxies is shown for the three redshift 
ranges. The distribution of the means with redshift is also shown. Reduction in the scatter of the parameter value with time is apparent. 
The fall in the mean value of asymmetry is more significant for the bulgeless disc sample. 



Figure 9. Asymmetry is plotted against half light radius for bulgeless disc and normal disc galaxies in the three redshift ranges. The 
solid line in each graph marks the linear relation followed by the full sample at that redshift range. Larger the radius, greater is seen to 
be the asymmetry value of the galaxy. Sources from the bulgeless disc sample, having on an average larger half light radii, are found to 
have higher asymmetry values. The typical error on the asymmetry value is ±0.08 and for the half light radius it is ±0.2 kpc. 


iLopez-Saniuan et alJ I 2 OIII : iBluck et al. I [mi). We do ob¬ 
serve that both bulgeless and normal disc galaxies have 
gained more than 50% of their total stellar mass in the past 
billion years. Although star formation within the galaxy 
can also lead to the increase in its stellar mass, its contri¬ 


bution is measured to be mu ch less compared to tha t from 
minor mergers and accret ion llOwnsworth et ^|2012| . l2014l : 
iMadau fc Dickinsonll20f3l . 

During the last Gyrs, galaxies are predicted to suf¬ 
fer a period of intense bombardment by minor satellites 
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Figure 10. Asymmetry is plotted against log of effective stellar mass density for bulgeless disc and normal disc galaxies in the three 
redshift ranges. The solid line in each graph marks the linear relation followed by the full sample at that redshift range. The two 
parameters are seen to be highly correlated for both morphological types. Higher the galaxy’s stellar mass density inside its effective 
radius, lower is the asymmetry of its stellar light distribution. The typical error on the asymmetry value is ±0.08 and for the log of 
effective stellar mass density it is ±0.05. 


dKhochfar & SillJ200f^: 

Hookins et alJ200fll: Feldmann et alJ 

l20inl: Oser et al.ll20inl: 

Ouilis & Truiillcj 2012|). These bom- 


bardments are expected to bring morphological changes 
in the disc population by building cl assical bulges and 
also giving rise to spheroidal galaxies dOes ch et a ll l2010l: 


Hookins et al.l2010l: 

van der Wei et al.ll2011 

: Cameron et al. 

2 OIII: IWeinzirl et al.ll201ll: iLaw et al.l 12012 

: iBuitraeo et al. 


2013fl . For the full sample studied, we report a considerable 


increase in the density of the inner region through the mea¬ 
sures of Sersic index, concentration as well as the effective 
stellar mass density for the past ~8 Gyrs. 


This increase is also observed in the fact that the frac¬ 
tion of disc galaxies with classical bulges is increasing from 
z ~ 0.9 to the present epoch. It increases from 60% at the 
highest redshift range (0.77-1.0) to 64% at the intermediate 
redshift range (0.4-0.77) and finally to 75% at the local red¬ 
shift range (0.02-0.05), indicating that some of the bulgeless 
discs are growing a classical bulge with time. Although in¬ 
ternal secular mechanisms driven by disc instabilities also 
lead to the increase in inner density, th ese methods can only 
lead t o the formation of pseudo-bulges jBinnev fc Tremaind 
19871: iKormendv &: Kennicuttlliool : iKormendv et al.ll2010l : 


Saha fc Gerhard 2013li . The evidence, thus, suggests that 


minor mergers and accretion are playing a significant role 
from 2 ~ 0.9 to the present epoch. 


There is substantial literature that establishes the de¬ 
pendence of galaxy properties on local environmental den¬ 
sity such that the higher the density of the local environ¬ 
ment, the more massiv e, dense, early - type and non-star- 
forming is the galaxy jDressleii Il984: iGomez et al.l l2003l : 


tormmg is tne galaxy iiuressieii i Vromez et ai.i i/uuj 

iBlanton fc Moustakasi[^009l : IScoville et al.ll2013l l. In the spe¬ 

cific case of discs, discs with cl assical bulges are rarely found 
in low density environments (IKormendv et aLll2010l l. Also, 
discs without classical bulges are expecte d to keep this way 
by b eing in more isolated environments dPeebles fc Nussen 
I 2 OIOII . By that argument, at 2 ~ 1, normal disc galaxies 
supporting a classical bulge and being denser, more massive 
and less star-forming are expected to be placed in denser 
environments as compared to the bulgeless disc galaxies. 

There is observational evidence that galaxies in denser 


environ ments show a more rapid increase of ga laxy size with 
redshift (ICooper et al.l[^012l : iLani et al.ll2013l l. Our results 
are in agreement, such that, normal disc galaxies show faster 
size evolution as compared to the bulgeless galaxies. Dense 
environments should also facilitate a stronger evolution in 
the inner region of galaxies as they are expected to undergo 
an in creased amount of accretion and interaction with satel¬ 
lites (IConseIicell2014l . and references therein). In our study 
also, in absolute terms, the increase in the inner density is 
observed to be more prominent in the case of the normal 
disc sample. Thus, environment appears to have strongly 
affected bulge growth over the past -^8 Gyrs. 

Overall, we have found that both internal and external 
mechanisms are involved in disc and bulge evolution. Exter¬ 
nal processes, in the form of minor mergers and accretion, 
appear to be playing a more effective role in growing classical 
bulges in relatively denser environments. 

Examining the inner region properties through both 
parametric and non-parametric measures provides consid¬ 
erable insight into the relative role of the processes involved 
in disc evolution. Further understanding can perhaps be ob¬ 
tained from studies in other wavelengths, especially the in¬ 
frared region, in which dust effects are minimised. 
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